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Observaciones operativas 


El autor presenta los resultados de un análisis re- 
gional de caudales máximos anuales realizado en 
tres áreas geográficas de Chile, que corresponden 
a su zona media o mediterránea. Los nombres, 
extensiones de cada área y número de estacio- 
nes hidrométricas procesadas son: (1) Región 
Metropolitana, con 15 403 km? y seis estaciones; 
(2) Región del Libertador Bernardo O'Higgins, 
con 16 387 km? y 16 estaciones, y (3) Región de 
Maule, con 30 296 km? y 28 estaciones. El autor 
señala que en esta zona la precipitación anual 
varía de los 250 a los 500 milímetros, y que sus 
sistemas hidrográficos tienen un régimen pluvial 
y pluvionival. 

Primera: todo análisis regional debe comenzar 
probando la calidad estadística de los registros 
hidrométricos por procesar del caudal máximo 
anual, pues de ello depende que los resultados de 
su análisis probabilístico sean confiables. Condi- 
ciones como aleatoriedad, independencia, homo- 
geneidad y estacionalidad deben ser verificadas 
estadísticamente, a través de diversas pruebas 
(Bobée $ Ashkar, 1991). El polemista observa que 
no se verificaron, mediante pruebas estadísticas 
generales o específicas, las condiciones citadas. 

Segunda: recientemente, Hosking y Wallis 
(1997a) han encontrado que los valores erróneos, 
los eventos dispersos (outliers), la tendencia, 
y los saltos o cambios en la media de los datos 
son reflejados en los momentos L de la muestra 
y, por ello, han sugerido una mezcla de los co- 
cientes L en un solo estadístico (D.), que mide la 
discordancia entre los cocientes L de cada registro 


y los promedios del grupo. Tal prueba permite 
encontrar registros que son discordantes con el 
grupo considerado como un todo. El polemista 
Opina que cuando menos la prueba de discordan- 
cias debió ser aplicada, para concluir respecto a 
la calidad estadística de los registros disponibles 
de caudal máximo anual. 

Tercera: el autor utiliza los tres registros más 
amplios de cada área y concluye respecto a su 
comportamiento que existe tendencia ascendente 
en todos y que, por lo tanto, los caudales de la 
zona mediterránea de Chile están mostrando 
valores más grandes en épocas recientes. En opi- 
nión del polemista, para que tal conclusión sea 
aceptable, primero se debió probar la homoge- 
neidad hidrológica de cada región, por ejemplo, 
a través de: (1) la versión corregida del Test de 
Langbein (Fill € Stedinger, 1995); (2) el Test H 
de Heterogeneidad (Hosking é: Wallis, 1997b), o 
bien, (3) la prueba basada en la dispersión de los 
coeficientes de variación (Wiltshire, 1986a, 1986b). 
Lo anterior, cuando se emplean los registros de 
caudal máximo anual para probar la similitud 
en la respuesta hidrológica de las cuencas que 
integran cada región. Sin embargo, Cunderlik y 
Burn (2002), y Ramírez et al. (2009) han sugerido 
que la fecha de ocurrencia media y la regularidad 
estacional de las crecientes pueden ser utilizadas 
como medidas de similitud, ya que implican 
una semejanza en características fisiográficas y 
meteorológicas, teniendo la ventaja de utilizar los 
registros hidrométricos para la verificación de la 
homogeneidad regional y no para su definición. 

Cuarta: el autor indica que las cuencas proce- 
sadas tienen un régimen pluvial y pluvionival, 
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es decir, que sus crecientes han sido originadas 
por dos mecanismos diferentes. En tales circuns- 
tancias, el polemista destaca que no se probó 
estadísticamente la presencia de dos poblacio- 
nes, para aceptar tal hipótesis o descartarla; por 
ejemplo, con base en pruebas no paramétricas, 
como las de Mann-Whitney, Kruskal-Wallis y 
Terry. 

Quinta: aceptando que no exista suficiente 
evidencia estadística para considerar los regis- 
tros procesados integrados por dos poblaciones, 
el polemista observa que no se ajustaron a los 
datos de caudal máximo anual los modelos 
probabilísticos que han sido sugeridos bajo 
precepto, como lo son las distribuciones: [1] 
Log-Pearson tipo III en EUA (WRC, 1977; 
Stedinger et al., 1993); [2] General de Valores 
Extremos (GVE) en Inglaterra en la década de 
1970 (NERC, 1975; Stedinger et al., 1993), y [3] 
recientemente la Logística Generalizada (Shaw 
et al., 2011). 

Sexta: el polemista también destaca que en 
la aplicación de los modelos probabilísticos 
utilizados no se indicó el método de ajuste que 
fue empleado, por ejemplo: momentos, máxima 
verosimilitud o momentos L. De mayor impor- 
tancia es la selección del mejor modelo proba- 
bilístico, lo cual actualmente se realiza a través 
del error estándar de ajuste o del error absoluto 
medio, que son medidas cuantitativas del ajuste, 
las cuales sí permiten discernir objetivamente 
entre un modelo y otro (Kite, 1977). 


Observaciones de concepto 


En el inciso titulado “Calidad de los ajustes”, el 
autor indica que en las tres regiones analizadas, 
la distribución de probabilidades Gumbel es la 
más adecuada y cita como referencia de soporte 
a Yue et al. (2006). En el final de tal inciso, el 
autor establece que se utilizó la distribución 
Gumbel debido a que es la más fácil de deter- 
minar y tiene menor número de parámetros por 
estimar. Posteriormente, el autor cita la referen- 
cia de Ónóz y Bayazit (2001) para indicar que 
tales autores encontraron que la distribución 
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Gumbel es la más recomendable al trabajar con 
caudales mínimos. 

Séptima: el polemista encuentra que la refe- 
rencia de Yue et al. (2006) trata la distribución 
de valores extremos bivariada, utilizando el 
modelo Gumbel como marginales; por lo tan- 
to, permite la estimación de la probabilidad 
conjunta de las variables caudal máximo y su 
respectivo volumen del hidrograma, así como 
del volumen, conjuntamente con la duración de 
la creciente. Por lo tanto, el autor está citando 
mal tal referencia, para referirse a la distribución 
Gumbel de dos parámetros de ajuste (ubicación 
y escala), la cual, por cierto, está incluida en el 
modelo GVE cuando su parámetro de forma (k) 
resulta cero. 

Octava: el polemista considera que el modelo 
GVE debió ser probado, a juzgar por los regis- 
tros del río Mapocho en Hacienda de Maipú, 
mostrado en la figura 3, y del río Claro en El 
Valle, expuesto en la figura 4, que son registros 
que probablemente tengan valores dispersos 
(outliers). 

Novena: el polemista encontró que la re- 
ferencia de Ónúz y Bayazit (2001) se refiere a 
series de duración parcial o caudales máximos 
superiores a un valor umbral, cuya ocurrencia 
de excedencias se modela con un proceso de 
Poisson, siendo la distribución exponencial la 
que define las magnitudes de tales caudales. La 
distribución exponencial es un caso del mode- 
lo Pearson tipo III y nada tiene que ver con la 
distribución Gumbel, excepto el contener en su 
fórmula al número e. Por lo anterior, el autor 
está citando mal tal referencia. 

El autor, al final del inciso denominado 
“Variabilidad regional”, indica que el caudal 
máximo anual en las tres regiones posee un 
aumento gradual en el tiempo, medido a través 
de la expresión del índice relativo de magnitud 
de crecidas (IRmoc). 

Décima: el polemista considera que por el 
hecho de no utilizar el caudal específico (m?/s/ 
km), es decir, por no tomar en cuenta el área 
de cuenca de cada estación de aforos, el [Rmc 
siempre definirá valores mayores que la unidad 
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al emplear los caudales máximos observados 
y también reportará magnitudes mayores al 
emplear predicciones, o caudales asociados a 
cierta probabilidad de excedencia, sobre todo 
en periodos de retorno superiores a la amplitud 
del registro. 

Décima primera: en opinión de polemista, la 
aplicación del Test U de Mann-Whitney por sub- 
tramos de registro es incorrecta, pues tal prueba 
busca las diferencias entre eventos o crecientes 
procedentes de cada mecanismo de formación o 
población diferente. La diferencia entre medias 
por tramos se busca con una prueba paramétrica 
t de Student o de Cramer, y no paramétrica de 
Wald-Wolfowitz. 

Décima segunda: el polemista pregunta, ¿por 
qué en la figura 4, a partir del año 2002, las 
gráficas de caudales máximos anuales siguen 
las líneas de tendencia lineal? 
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RESPUESTA 


En primer lugar, es importante destacar la 
prolijidad de los alcances que hace el polemista 
y que se haya interesado en la temática que ex- 
pusimos con mis colegas en el paper en cuestión. 
En este marco es también importante señalar 
que el único objetivo que tuvo la investigación 
hecha fue intentar determinar si es que se había 
verificado algún cambio en el patrón de com- 
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AL POLEMISTA 


portamiento de las series de caudales máximos 
en diversas cuencas de la zona central del país. 
Para ello se utilizaron algunas funciones de 
distribución de probabilidad (FDP), escogién- 
dose la de Gumbel por tener una buena calidad 
de ajuste y estimando a partir de la misma los 
valores de caudal para los periodos de retorno 
de 10, 20, 30, 40 y 50 años. A partir de esto, fue 
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posible señalar que ha existido un aumento de 


los caudales máximos en los últimos periodos, 


lo que es refrendado incluso por medidas tan 


básicas como lo es la comparación gráfica y las 


tendencias de los datos. Por tanto, no hay mayor 


análisis sobre lo expuesto que lo ya dicho. 


A continuación se da respuesta a los alcances 


del polemista: 
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El boletín 17B del WRC señala que los 
caudales máximos anuales pueden ser 
tratados como eventos independientes, 
aunque las pruebas para verificar este 
supuesto indiquen lo contrario. En cuanto 
a la calidad de los datos empleados en el 
análisis, la Dirección General de Aguas de 
Chile realiza un exhaustivo análisis de los 
datos que pone a disposición de los usuarios, 
los cuales son validados de 2 a 3 años 
después de la ocurrencia de los mismos, 
con lo cual se verifica el comportamiento 
de las curvas de descarga. Esto define que 
los datos son de calidad y se ajustan a la 
realidad descrita en cada cuenca. Además, 
dado que el largo de los registros es superior 
a 10 años, los datos se consideran aceptables 
para realizar un ajuste de FDP. 

El ajuste de las FDP no consideró un análisis 
con base en los L momentos. Además, el 
análisis se realizó por estación, dadas las 
características particulares que posee el 
territorio chileno, que hace muy singulares 
a cada una de sus cuencas (caen más de 
5 000 m en menos de 100 km de recorrido; 
algunas poseen régimen pluvionival y otras 
sólo de tipo pluvial; las pendientes medias 
del cauce principal son muy distintas, 
dependiendo de si la cuenca se ubica en 
la cordillera de los Andes o de la costa, 
etc.). Por ende, cualquier agrupación de 
los datos en términos regionales estaría 
generando la necesidad de establecer 
análisis muy distintos a los efectuados, 
por la multivariabilidad que habría que 
considerar. En función de lo anterior, la 
aplicación del índice Di no tendría sentido, 
ya que éste se utiliza para la identificación 


de medidas discordantes dentro de un 
análisis grupal, que para este caso no 
corresponden. 

Sólo se seleccionaron los registros más 
amplios de cada área, pues existe mucha 
variabilidad en la longitud de las series, 
y se asumió que ellas podían representar 
al conjunto de una región administrativa, 
no hidrológica ni estadística. El polemista 
insiste en entender, al parecer, que las 
regiones son de tipo hidrológico y eso no 
corresponde al artículo. 

No se planteó la hipótesis de que los 
caudales estudiados proviniesen de dos 
regímenes distintos, dado que lo que se 
desea probar es si ha habido cambios 
temporales en el comportamiento de los 
caudales máximos para cada cuenca en 
estudio y de forma particular; esto es, si 
ellos han aumentado en los últimos años. 
Lo expuesto por el polemista tendría 
mayor asidero (clasificando en cuencas 
pluviales y pluvionivales), si este estudio 
correspondiese a una descripción de la 
zona de estudio en términos de análisis 
regionales. Asimismo, todas las cuencas 
en estudio corresponden a pluviales y 
pluvionivales (Cortés et al., 2011; Bravo éz 
Garreud, 2011). 

Los modelos probados en el estudio son 
utilizados comúnmente en Chile. Además, el 
estudio no pretendía encontrar la mejor FDP 
para las estaciones, y dado que los valores 
del coeficiente de determinación y prueba 
de Kolmogorov-Smirnov fueron buenos, no 
fue necesario ajustar otro modelo de FDP 
conforme a los objetivos del estudio. 

El ajuste de los modelos se realizó a través 
del método de los momentos. En cuanto a 
la selección de los modelos, los indicadores 
que señala el polemista presentan el error en 
la misma unidad que la variable, por lo que 
su interpretación es sencilla, pero no son 
objetivos. Asimismo, Chai y Draxler (2014) 
señalan que los indicadores (RMSE Y MAE) 
sólo enfatizan un aspecto del error, por lo 
que lo ideal es utilizar una combinación 
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de pruebas de bondad de ajuste para la 
selección del modelo. 

7. En términos estrictos, el polemista tiene 
razón; sin embargo, lo que se quiere destacar 
es la ductilidad y aplicación de la FDP de 
Gumbel y no necesariamente la aplicación 
específica que de la misma hicieron los 
autores Yue et al. (2006). 

8. Los modelos probados en la investigación 
tuvieron un buen ajuste con base en el 
coeficiente de determinación (95% para 
Gumbel y 98% log-normal) y la prueba 
de Kolmogorov-Smirnov no rechaza la 
hipótesis nula en ambos modelos, por lo 
que se acepta que el ajuste es bueno. El 
ajustar otro modelo sólo sería recomendable 
en el caso de que los modelos probados no 
tuvieran un ajuste aceptable. 

9. Existe coincidencia con el polemista. 

10. El polemista tiene razón, porque es 
justamente lo que se pretendía hacer; esto es, 
definir el valor uno como el valor referencial 
para establecer los análisis comparativos en 
el marco del índice propuesto. 

11. La U de Mann-Whitney es una instancia 
que permite hacer comparaciones en datos 
no paramétricos y a partir de las medianas. 
El objetivo de su aplicación fue determinar 


si había diferencias entre las series 
consideradas y si esto permitía corroborar 
la hipótesis de que los caudales punta se 
habían incrementado, lo cual se consiguió 
no sólo con esta prueba, sino que también 
con el análisis de tendencia y el gráfico. De 
igual manera, las pruebas sugeridas por el 
polemista presentan problemas similares a 
los criticados a la U de Mann-Whitney. 

12. Corresponde a un error gráfico que sólo se 
manifiesta en la región de O'Higgins. 
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